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Abstract. Electromagnetic streaming instabilities of multicomponent collisional magne- 
tized accretion disks are studied. Sufficiently ionized regions of the disk are explored where 
there is strong collisional coupling of neutral atoms with both ions and dust grains simulta- 
neously. The steady state is investigated in detail and the azimuthal and radial background 
velocities of species are calculated. The azimuthal velocity of ions, dust grains, and neutrals 
is found to be less than the Keplerian velocity. The radial velocity of neutrals and dust 
grains is shown to be directed inward of the disk. The general solution for the perturbed 
velocities of species taking into account collisions and thermal pressure is obtained. The 
effect on the collisional frequencies, due to density perturbations of charged species and neu- 
trals, is included. It is shown that dust grains can be involved in the fast electromagnetic 
perturbations induced by the ions and electrons through the strong collisions of these grains 
with neutrals that in turn have a strong collisional coupling with the ions. The dispersion 
relation for the vertical perturbations is derived and its unstable solutions due to different 
background velocities of ions and electrons are found. The growth rates of the streaming 
instabilities considered can be much larger than the Keplerian frequency. 
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1. Introduction 

Protostellar and protoplanetary accretion disks are known to contain submicron- and 
micron-sized solid particles (dust grains) (e.g., Beckwith & Sargent 1996; Isella et al. 2006; 
Besla & Wu 2007; Quanz et al. 2007; Pinte et al. 2007). The presence of dust grains 
is common in these and other astrophysical and cosmic objects, such as, molecular and 
interstellar clouds, planetary and stellar atmospheres, planetary rings, and cometary tails 
(Spitzer, Jr. 1978; Goertz 1989; Havnes et al. 1996; Wardle & Ng 1999; Rotundi et al. 
2000). Dust is contained in various plasmas of both research and technological importance 
(Merlino & Goree 2004). In space, dust grains vary widely in sizes, from submicrons up to 
1 cm and more. Dust particles in both space and the laboratory acquire an electric charge 
due to cosmic-ray, radioactive or thermal ionization of the ambient gas, photoionization, 
absorption of electrons and ions from the background plasma, the presence of electron and 
ion currents, thermoionic emission, and secondary ionization (Hayashi 1981; Whipple 1981; 
Meyer- Vernet 1982; Chow et al. 1993; Barkan et al. 1994; Melzer et al. 1994; Mendis & 
Rosenberg 1994; Thomas et al. 1994; Horanyi 1996; Fisher & Duerbeck 1998; Fortov et al. 
1998). The polarity of charged grains in space and the laboratory can be negative or positive 
(Meyer- Vernet 1982; Chow et al. 1993; Mendis & Rosenberg 1994; Horanyi 1996; Horanyi 
& Goertz 1990; Havnes et al. 2001; Ellis & Neff 1991). The magnitude of grain charges in 
the laboratory can be very large (up to a few orders of magnitude of the electron charge). 
The polarity of charged grains depends on their size; large grains arc usually negative and 
small grains positive (e.g., Chow et al. 1993; Mendis & Rosenberg 1994). 

The dynamics of dust grains in protostellar and protoplanetary accretion disks deter- 
mines the physical processes leading to the formation of planets. It is generally accepted that 
small solids (1 cm - 1 m) can be formed due to chaotic motion and sticking of submicron- 
and micron-sized dust grains strongly collision coupled to the surrounding gas. These small 
solids in turn can form the kilometer-sized bodies (planetesimals) , believed to be the building 
blocks for planet formation. 

Various physical mechanisms have been discussed for the formation of planetesimals, 
such as, coUisional coagulation (Adachi et al. 1976; Cuzzi et al. 1993; Weidenschilling 1995), 
gravitational instabihty (Safronov 1969; Goldreich & Ward 1973; Yamoto & Sekiya 2004), 
and the trapping of particles by large-scale persistent vortices (Barge & Sommeria 1995; 
Lovelace et al. 1999; Klahr & Bodenheimer 2003; Johansen et al. 2004; Petersen et al. 
2007). 

Recently the streaming instability arising from the difference between velocities of dust 
grains and gas has been suggested as a possible source contributing to planetesimal formation 
(Youdin & Goodman 2005). These authors have considered the hydrodynamic instability of 
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the interpenetrating streams of dust grains and gas coupling via drag forces in the Keplerian 
disk. The particle density perturbations generated by this instability could seed planetesimal 
formation without self-gravity. 

In the papers cited above, which deal with the dynamics of dust grains in accretion disks 
in connection with the problem of dust coagulation and planet formation, hydrodynamic 
models were applied in which the dust grains were considered as a neutral fluid or ensemble 
of individual particles interacting with the surrounding neutral gas. Thus, electrodynamic 
processes were not involved. However, dust grains are generally charged (sec references 
cited above). In the surface layers and near the star, the protostcllar and protoplanctary 
accretion disks are multicomponent, containing electrons, ions, charged grains, and neutral 
gas. In the sufficiently dense inner regions, the disk matter likely contains only charged 
dust grains and gas (e.g., Wardle & Ng 1999). Dust grains can be of various forms. It is 
clear that electromagnetic forces could play an important role in dust coagulation and in 
collective processes within the disk involving the dust grains. The electromagnetic turbulence 
favors the coalescence of charged dust grains. The collisional coagulation due to electrostatic 
interaction of dust grains can lead to the formation of fractal aggregate structures (Matthews 
et al. 2007). These structures could build larger-sized solid particles in the process of planet 
formation. 

The necessity of taking into account electromagnetic phenomena in the disk dynamics 
is confirmed by observations, which show that astrophysical disks arc turbulent (Hartmann 
2000; Carr ct al. 2004; Hcrsant et al. 2005). It is also known that accretion disks arc threaded 
by the magnetic field (Hutawarakorn & Cohen 1999 and 2005; Donati et al. 2005). In the 
latter, a direct detection of the magnetic field in the protostellar accretion disk FU Orionis, 
including the innermost and densest parts of the disk, is reported. The surface magnetic 
field was observed to reach strengths of about 1 kG close to the center of the disk and 
several hundred gauss in its innermost regions. The magnetic field is grossly axisymmetric 
and has both vertical and significant azimuthal components. The results of Donati et al. 
also suggest that magnetic fields in accretion disks could trigger turbulent instabilities that 
would produce enhanced radial accretion and drifts of ionized plasma through transverse 
field lines. 

The electromagnetic dynamics of accretion disks is determined by physical parameters, 
such as, the magnetic field, density, temperature, the ionization degree, composition, and 
mass and size of charged dust grains, all of which vary significantly through the whole 
disk. So these would determine the disk dynamics on a local scale. The protostellar and 
protoplanctary accretion disks as well as cold molecular and interstellar clouds are weakly 
ionized (Wardle & Ng 1999; Norman & Heyvaerts 1985; Jin 1996; Bergin et al. 1999; Balbus 
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& Terquem 2001; Desch 2004; Tscharnuter & Gail 2007). Dust grains are, possibly, the 
primary charge carriers in very dense, cool nebulae (Wardle & Ng 1999; Blaes & Balbus 
1994). 

One of the sources of electromagnetic turbulence in accretion disks may be electromag- 
netic streaming instabilities emerging due to the different velocities of species in equilibrium. 
We have studied this issue (Nekrasov 2007) for cold magnetized disks. The dust grains 
were considered to be magnetized and strong coUisional coupling of neutrals with one of 
the charged species could be included. The thermal pressure was not taken into account. 
New compressible instabilities were found, with growth rates much larger than the Keplerian 
frequency. In a subsequent paper (Nekrasov 2008), the general theory for electromagnetic 
streaming instabilities in multicomponent weakly ionized regions of accretion disks embedded 
in the magnetic field has been developed. The compressibility, anisotropic thermal pressure, 
and collisions of charged species with neutrals were taken into account. The equilibrium state 
was found for the case in which the ions do not influence the motion of neutrals. However, 
the neutrals are assumed to have weak or strong collisional coupling with dust grains. In 
the perturbed state, the neutrals have been considered as immobile when the perturbation 
frequency is much larger than the collisional frequency of neutrals with charged species. A 
dispersion relation has been derived in the most general form for an arbitrary direction of 
the perturbation wave vector and for arbitrary strength of thermal effects. In particular, 
we found solutions of the dispersion relation for a case appropriate to definite regions of the 
protostellar and protoplanetary disks, where the electrons can be considered as thermal and 
magnetized, and ions and dust grains as cold and unmagnetized. These solutions describe 
new instabilities of the weakly ionized disks due to collisions and differences between the 
stationary velocities of different charged species. 

In the present paper, we study the electromagnetic streaming instabilities in the suf- 
ficiently ionized regions of the protostellar and protoplanetary magnetized accretion disks, 
where the neutrals have strong collisional coupling simultaneously with ions and dust grains 
both in the equilibrium and in the perturbed state. Dust grains and light charged com- 
ponents are respectively unmagnetized and magnetized, i.e., their cyclotron frequencies are 
respectively much smaller or much larger than their orbiting frequencies. The dust grains are 
treated as monosized with a constant charge. For the coUisional regime under consideration, 
we find the azimuthal and radial stationary velocities of species. We consider the effect of the 
change in collisional frequencies due to density perturbations of species. Taking into account 
anisotropic thermal pressure, we derive the general expression for the perturbed velocity of 
any species that contains also the perturbed velocity of other species due to collisions. In 
the cold limit for horizontally elongated perturbations, we obtain the general solutions for 
the perturbed velocities of ions and dust grains incorporating their mutual influence on each 
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other via collisions with neutrals. Further, we consider fast processes in which the electro- 
magnetic dynamics of dust grains is determined by collisions with neutrals but not by their 
own motional frequencies. Finally, we investigate unstable perturbations due to different 
azimuthal velocities of ions and electrons. 

The paper is organized as follows. In § 2 the basic equations are given. In § 3 we 

discuss in detail the coUisional equilibrium state. The general solutions for the perturbed 
velocities and densities of species including thermal effects are obtained in § 4. The perturbed 
velocities of species for vertical perturbations in the cold limit are considered in § 5. Solutions 
for the perturbed velocities of species when the dynamics of dust grains is determined by 
collisions with neutrals are given in § 6. In § 7 we calculate the perturbed electric current. 
The dispersion relation in the general case is derived in § 8. The dispersion relation in the 
case of strong coUisional coupling of dust grains with ions through collisions with neutrals is 
considered in § 9. Its unstable solutions are found in § 10. Discussion of the obtained results 
is given in § 11. The main points of the paper are summarized in § 12. 



2. Fundamentals 

The fundamental equations in the inertial (nonrotating) reference frame are the follow- 
ing: 

^ + V,- • Vv,- = -VU - ^-p- + -^E+-^v,- X B-i^jn (v,- - v„) , (1) 



dt 

the momentum equations, 
the continuity equation. 



V-P, 


+ E+- 




rrij r 


-VU- 









^ + V„ • Vv„ = -VU - l^nj (Vn - Vj) , (2) 



_^ + V-n,- „v,- „ = 0, (3) 



VxE^-lf. (4) 

V X B =— j, (5) 

c 

and Maxwell's equations, where j = qjUjYj. Here the index j — e,i,d denotes the elec- 
trons, ions, and dust grains, respectively, and the index n denotes the neutrals. In Eqs. 
(l)-(5) Qj and mj^n are the charge and mass of species j and neutrals, Vj^n is the hydro- 
dynamic velocity, rij^n is the number density, Pj and P„ are the thermal pressure tensor 
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of charged species j and the thermal pressure of neutrals, respectively, Ujn — ^jrjnnnn 
i^nj = ^jniTT-jnj) is the coUisional frequency of charged species (neutrals) with neutrals 
(charged species), where 7j„ =< av >jn /{jnj + m„) (< av >j„is the rate coefficient for 
momentum transfer), U = —GM/R is the gravitational potential of the central object hav- 
ing mass M, R = (r^ + z'^Y^'^, G is the gravitational constant, E and B are the electric 
and magnetic fields, and c is the speed of light in vacuum. We consider wave processes with 
typical time scales much larger than the time the light spends to cover the wavelength of 
perturbations. In this case one can neglect the displacement current in Eq. (5), what results 
in quasineutrality both in electromagnetic and purely electrostatic perturbations. We use 
the cylindrical coordinate system {r,9,z), where r is the distance from the symmetry axis 
and z is the coordinate from the midplane of the disk. The magnetic field B includes the 
external magnetic field ^oext of the central object and/or interstellar medium, the magnetic 
field Bocur of the background current, and the perturbed magnetic field. The pressure tensor 
Pj has the form Pj=p±jl+{p\\j —p±j)hh, where p±j and p\\j arc the pressure across (_L) and 
along (II) the magnetic field B, I is the unit tensor, and b is the unit vector along B. 

Equations (1) and (2) take into account only the collisions between charged species and 
neutrals. In the weakly ionized objects these collisions are dominant. The system of Eqs. 
(l)-(5) can be applied in the coUisionless regime, if we put = u^j = 0, as well as in the 
regimes of weak and strong coUisional coupling of neutrals with charged components. In Eqs. 
(l)-(3) we do not take into account ionization and recombination processes. We comment 
on this point in § 11. Self-gravity is not included in the present paper. 



We suppose that Boe^t is axisymmetric and -Boeeit = 0. Such a configuration is typical 
of the magnetic field of the central star and/or of disks threaded by the vertical interstellar 
magnetic field. We neglect the radial component of the background magnetic field, consid- 
ering the regions of the disk, where vertical components of both the external magnetic field 
and the magnetic field induced by the stationary azimuthal current are dominant. 

Let us consider axisymmetric stationary flows of species j. We suppose that the vertical 
stationary velocity Vjo^ is equal to zero. Then the r and 9 components of Eq. (1) in the 
equilibrium take the form 



3. Equilibrium 




vl 



r 



+ CO, 



(6) 
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Here cucj — qjBo/mjC is the cyclotron frequency of species j {Bq — Bo^ext + -Bozcur), is 
the coUisional frequency in the equihbrium, and 

dU 1 dp±_Qj 

Fjor = —rrij-^ ^ \- QjEor, 

or rijo or 

where Eor is the background electric field. 

Equations (6) and (7) determine the stationary velocity v^o of different charged compo- 
nents due to the action of the electric (radial) , magnetic (vertical) , and gravitational fields as 
well as the thermal pressure and coUisions with neutrals. Due to the latter effect the radial 
velocity emerges [Eq. (7)]. Solving Eqs. (6) and (7), we will consider two cases in which there 
are magnetized, cUcjVjoe ^ "^jorei/^' unmagnetized, uicjVjoQ -C v'jQrd/''^^ charged species. 
The condition Ucj 3> (<^)fjoe/'" denotes that the cyclotron frequency is much larger (smaller) 
than the orbital frequency [we use below, as usual, the term "magnetized" ("unmagnetized") 
also in the case in which cUcj ^ {^)^jn]- 

For magnetized species we have the following solutions of Eqs. (6) and (7): 

^cj FjOr ^jn ^cj^jn /q\ 

^cj + ^jn <^cj + ^jn ^cj + ^jn 

^jOr - 2 , ,,2 ^ + 2 , ,,2 ^«0e + 9 , ,,2 ^nOr- [^) 

Here, for simplicity, the index at has been omitted. 

We will consider the case in which simultaneous coUisions of neutrals with two charged 
components are important. We suppose that one of these components is the heavy unmagne- 
tized dust grains d and another component, which will be denoted by the index j, represents 
the ions or the light magnetized dust grains (the latter case is appropriate for systems con- 
taining only two dust components). Due to the small electron mass the collisions of neutrals 
with electrons are neglected. Below in this section, the corresponding condition is given. 
Then, the equation of motion (2) for neutrals takes the form 

= VnA'^nQr - VjOr) - l^nd{VnOr " VdOr), (10) 

Here i'n3,d — ^nj,d- unmagnetized dust grains Eqs. (6) and (7) have the form 

d^dor '^m _ FdOr , X 

— T^dn[VdOr — VnOr), U^j 



2dr r rrid 
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1 drv. 



doe 



r dr 



VdOr — —^dn{vdae — Vnoe)- 



(13) 



To solve the system (10)-(13) we adopt, for simplicity, that Vn,doe ~ and Vn,dOr ~ 
r~"2 Then we can substitute the derivatives in Eqs. (10)-(13) by the algebraic terms. We 
find further the expressions for w„or ~ "^jOr and t;„oe ~ Vj^g, using Eqs. (8) and (9), and 
substitute them in Eqs. (10) and (11). From Eqs. (11) and (13) we can find the radial 
velocities VnOr and Vdor- 



VnOr 



jOr 



+ UJcj I - 



[1-x), 



VdOr = -1^1^(1 - -), 



X 



(14) 
(15) 



where x — Vdoehnoe and the index * denotes the modified collisional frequencies I'^j^a 
^nj,d/{^ — CKi) and = T^dn/i^ — Oil)- The parameter aj in Eq. (14) is 



1 + 6,- = 1 + 



(16) 



After substitution of expressions (14) and (15) into Eqs. (10) and (12) we obtain two 
equations containing only the azimuthal velocities of neutrals and dust grains (except the 
terms proportional to 0:2): 



■'noe 



0L1 



nOr 



+ 



"nOr 



+ 



7/* f -2 



nd 



+ UndT (1 - X) 



V. 



dn 
X 



+ 



7/* I ? 

nj cj 



nd 



^doe 



Ol2 



+ 



dOr- 



l^dnKj^cj ( FjOr 



(17) 



uJcPnoe]-i^dnr{l-x)[^+^-^]. 



^nd 



X 



(18) 



For simplicity, we have retained in these equations the radial velocities in terms proportional 
to a2. Multiplying Eq. (17) by x'^ and equating it to Eq. (18), we obtain the following 
equation: 



x 



2 FnOr FdQr 



v., 



nj 



rrir. 



jOr 



2 I / 2__n£_cj_ 



x'^Kd + ^dn 



ajVnoe, 



+r (1 — x) 



X 



7/* , i2 



nd 



C3 



jn 



.a,^+x''a2^ = 0. 
r r 



(19) 
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In the paper (Nekrasov 2008), we have considered the case of the weak coUisional cou- 
phng of neutrals with ions in the steady state, which is appropriate for weakly ionized disks. 
The parameter bj [sec Eq. (16)] for ions has been supposed to be much less than unity, in 
which case ^ 1. This condition can be satisfied for magnetized, ujd ^ z^m, or unmagne- 
tized, Uci ^ i^in, ions, if the orbital (Keplerian) frequency of neutrals Vnoe/f^ is less than the 
coUisional frequency of neutrals with ions Uni- In the case uJci ~ ^in it should be z/„j -C Vnoe/r. 
In the present paper, we will consider the opposite case in which 6j ^ 1, i.e., 

» (20) 

^cj ' jn ' 

Thus, we explore here the regions of the disk, where the medium is sufficiently ionized and 
where a strong coUisional coupling between neutrals and species j (ions or light dust grains) 
takes place. 

We can find the analytical solution of Eq. (19) in the case in which the azimuthal 
velocities of neutrals and dust grains are close to one another, i.e., 

a; = 1 + 5,5 < 1. 



This condition is easily satisfied in disks, where the coUisional coupling of dust grains and 
neutrals is strong (see below) . We neglect further the contribution of terms proportional to 
S in the first and third terms on the left-hand side of Eq. (19) as well as the contribution 
of two last terms in this equation. The corresponding conditions, taking into account the 
inequality (20), have the form 



2 2 
r 



e 2 -L / 

01' > a2-j 



nOr 



;2a; 



?;2 



). 



(21) 



where ~ I'dni^dn -\- ^nd)- We suppose that hjVdn ^ ^nd- This condition is, obviously, 
satisfied. Conditions (21) imply strong dust-neutrals coUisional coupling. Then we obtain 
from Eq. (19) the following solution for 5: 



"nOr 



rrid 



H i^ni + ^nd + ^dn] 



+ ^^cjVnOe 



(22) 



The second equation connecting S and Vnoe is obtained from Eq. (18) under conditions 
(20) and 

(23) 



S^dn^dn » ^ + «2 



■^dOr 
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This inequality also implies strong coUisional coupling of dust grains with neutrals. Then 
we have 



^f^dn'^dn = \ + ^cjVnOe ■ (24) 

rrid Ujn \mj ' J 

We can find fi:om Eqs. (22) and (24) the expressions for Vnoe and S. The velocity Vnoe 
is equal to 

FjOr ^jn I FnOr ^nd FdQr 



- i^cjVnoe = + — + • (25 

rrij Unj V l^dn rUd 



It can be shown that the condition 



must be satisfied for the solution (25) to have the present form. It is obvious that this 
inequality can easily be realized, if bjVdn ^ ^njd- The value 5 we find from Eqs. (24) and 
(25), 

1 ( FnOr ^nj ~l~ ^nd FdOr 



^^^dn = + — • (26) 

Unj V l^dn rUd J 

The azimuthal velocity of dust grains is equal to 

Vdoe = Vnoe{^ + S)- (27) 

The radial velocities of neutrals and dust grains are found from Eqs. (12) [or (14)] and (15): 

FdOr . * 
VnOr = Sru^n: 

^dn-nid 

Vdor = -Sri^dn- (28) 

Let us now consider the velocities Wj06»,r- Substituting the velocities Vnoe,r found above 
into Eqs. (8) and (9), we obtain Vjoe ~ Vnoe and Vjor ^ VjQg,Vn,dor- 

The electron velocities under condition lu'^ :$> v^^ are the following [see Eqs. (8) and 

(9)]: 

I'eOe = — 1 '^VnW —VnOri 



FeOr 








1 :rVnoe 






ce 




^en F(,Q.f 








1 Vnoe ■ 








ce 



VeOr = 1 VnOe + ^^^nOr- (29) 
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We do not take into account the collisions of neutrals with electrons. It can be shown 
that the corresponding condition has the form 



, ,2 I ,,2 • 



We suppose that this inequality is satisfied. 



4. Perturbed velocities and densities of species: General solutions 

In the present paper, we do not treat perturbations connected with the background 
pressure gradients. Thus, we exclude from our consideration the drift and internal gravity 
waves. We take into account the induced reaction of neutrals on the perturbed motion of 
charged species. This effect is important, if the ionization degree of medium is sufficiently 
high. We also include the effect of perturbation of the coUisional frequencies due to density 
perturbations of charged species and neutrals. This effect emerges when there are different 
background velocities of species, as occurs in accretion disks. Then the momentum equations 
(1) and (2) in the linear approximation take the form 

^+v,o-Vv,i+v,rVv,o = — ^Vn,i + ^Ei + ^v,oxBi + ^v,ixBo+C,i, (30) 
at 'mj'rijo rrij rrijC rrijC 

+ Vno • Vv„i + V„i • Vv„o = -%n hC„i, (31) 

In these equations, the coUisional terms C^i and C„i are equal 

Cjl = -I'jn^jl + Rjl, 

where = ■jjn'nT'nnno, = J2j^njy ^nj = JnjiTT'jn'jo- The index 1 denotes the quantities 
of the first order of magnitude. Here Pqj = 7±jP±ojI + {l\\jP\\oj — l±jP±oj)^^, l±j = 2, 
7||j = 3 for magnetized particles ujcj 3> z/j„. For charged species wc use here the double 
adiabatic dependence of temperature perturbations on the density perturbations analogous 
to the double adiabatic law of Chew et al. (1956) for the coUisionless plasma. For un- 
magnetized charged particles, cucj <S I'jn, we have Poj = JjPoj^, Ij — 5/3. In the case of 
Boltzmann's distribution of the number density due to movement of particles along the 
background magnetic field we have — 1. The neutrals are also considered as adiabatic. 
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Then V Pni/mnUnQ — i'TnVn„i/n„o, where v'^^ — ^nT^vaj'^ni Tno is the unperturbed temper- 
ature of neutrals. 



The term Rji, 



Vnl H (VnO - Vjo) 



describes the back reaction of charged components due to influence of the perturbed motion 
of neutrals and the term R„i, 



3 



Vji H (Vjo - v„o) 



leads to the induced perturbation of the velocity of neutrals. The coUisional terms propor- 
tional to n„i and nj\ take into account the dependence of the coUisional frequencies Vj^ and 
on the density perturbations of species. 

In this section, we will consider the general case in which the wave vector k = (/c^, m^, 
of perturbations has three components. We provide the normal mode analysis, assume 
the perturbations to be proportional to exp(iA;,.r + imQ + ik^z — iujt), and use the local 
approximation |k|r ^ 1. Then we find the Fourier amplitudes of the density perturbations 
from the linearized version of the continuity equation (3) 

f^j,nlk kVj 

where ujj^n = '-^ ~^'^j,no, k = {k,a;}. Substituting n„i into Eq. (31) and solving this equation 
in the Fourier representation, we will find the velocity v^i^, whose components and the 
expression for Unik are given in the Appendix A. 

In the protostellar and protoplanetary disks, the collisions of neutrals with charged 
species are frequent enough that the condition uj,,n ^ is satisfied (for systems with the 
weak collisions it should be supposed that uin ^ ^n)- Then, using the equation for Unik 
from the Appendix A and the expression for R„i, we obtain 

Dn = ^UJ^n V C-^- 32 

rino ^ ■' rijo 

In the present paper, we will further neglect the thermal effects of species when con- 
sidering some specific cases. For neutrals, the necessary condition for this cold limit has 
the form 6imn ^ 1, where l,m = r,9,z. In this case, the perturbed velocity and density of 
neutrals given in the Appendix A and Eq. (32), correspondingly, will be equal to 

Vnlfe = R-nlfc (33) 
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and 



rino 



I 



E 



V, 



^jlfc 



(34) 



Let us substitute the expressions (33) and (34) into Rjifc. Then the coUisional term Cjifc can 
be written in the form 

Cjife = —i^j^jik + Qjik, (35) 

where 



vufe H (Vio - Vjo) 



Here the index / denotes the charged components. The coUisional term Q^i describes the 
reaction of the species j on the coUisions of neutrals with another charged components. 

In the weakly ionized plasma the frequency of collisions of neutrals with charged com- 
ponents can be small compared to the frequency of perturbations, which depends on the 
disk parameters, a;„ ^ z/°. In this case, the induced dynamics of neutrals can be neglected, 
i.e., Vj PS 1/°^ and Q^i^ ^ 0. Such a situation was considered by Nekrasov (2008). In the 
sufficiently ionized plasma the opposite case can be realized in which i/° ^ Then the 
coUisional frequency Uj takes the form 



7/0 




(36) 



We see that when there is only one charged component with which the neutrals collide, the 
inertia of neutrals is added to the inertia of this charged component, i.e., Pj = —iujni^j^/ ^nj)- 
At the same time the value Q^i = 0. Thus, in other respects the plasma stays formally 
coUisionless. However, the presence of other charged components with which neutrals collide 
with the sufficiently high frequency, v^i 3> cu^, changes the situation dramatically: The 
plasma stays coUisional and the value Q^i contributes to the momentum equation of species 
j. Such conditions are realized in the ionized regions of protostellar and protoplanetary 
disks. This case will be considered in the present paper. 

The solution of Eq. (30) with Cjik defined by Eq. (35) and the expression for rijik are 
given in the general form in the Appendix B. 



We will also neglect the thermal effects for the charged species. The necessary condition 
is <S 1, where Z, m = r,9,z. We consider further the case in which the neutrals collide 
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simultaneously with ions and dust grains. As a result, the mutual influence of perturbations 
of ions and dust grains via perturbations of neutrals arises. The corresponding system of 
equations describing this process is given in the Appendix C. 

5. Perturbed velocities of species in the case ^ = kg = 0: General 



We have obtained the exact solutions of the system of equations for vnrft and Vdir,e given 
in the Appendix C for the horizontally elongated perturbations when kz ^ kr — ke — 
These solutions and the general solutions for Vi^^iz a-re given in the Appendix D. 

If we set //j = //d = in Eqs. (C.l) and (C.2), i.e., if we neglect the coUisions of neutrals 
with ions and dust grains (a;„ > then we obtain A,d^^i,dir = -Hi^dir, Di^dVi,die = -Hi^die, 
and Vi^diz = ii/^vi,d)Gi,diz- In this case, the ions and dust grains move independently from 
each other under the action of the electromagnetic forces. However, in the case u;„ <^ the 
movement of dust grains depends on the movement of ions, and vice versa (see below). 

The general solutions for the perturbed electron velocity in the case kz ^ 0, kr = kg = 
are given in the Appendix E. One can see that the collisions of neutrals with ions and dust 
grains influence the collisions of electrons with neutrals. 

6. Solutions for v^i^ in the case of neglect the own dynamics of dust grains 

The own frequencies of dust grains, the plasma and cyclotron frequencies, are very low 
because of their large mass. Therefore, at the absence of collisions of species with dust grains 
the latter can stay immobile in the electromagnetic perturbations having sufficiently high 
frequency. However, for example, through collisions with neutrals, the dust grains can be 
involved in the fast electromagnetic perturbations induced by the ions and electrons when 
the ions have also strong collisional coupling with neutrals. 

Let us neglect in Eq. (D.3) the terms proportional to cudi and 00^2 under conditions 



Note, that for frequencies of perturbations much less than the collisional frequencies, Ui^d S> 
cui^d, these conditions can be written in the form i/^ ^ (1 + i'i/oJci)uJdi,2- We have taken into 
account that Ud ^ 0,^. Then the value A will be equal to A — — A, where A = 
In the case (37), one can also neglect the action of the electromagnetic forces on the dust 



solutions 




(37) 
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grains. Then the expressions (D.l) and (D.2) take the form, 



Viie — {—^ud^ciGiir + ^ujGiie) ■ (38) 

For the dust grains we have 

l^i U C 

Vdw = — —Viw - -^k^weGiu. (39) 
The vertical velocities vnz and Vdiz are equal in the case under consideration, 

Kvdu = -l^iGiu- (40) 

We can find now the perturbed electron velocity. Substituting the velocities Vi^di defined 
by Eqs. (38) and (39) in the expressions for the components of the electron velocity given 
in the Appendix E, we obtain 

DeVelr — —Heir + ^ {^ui^ue ~ i^ud^ce^ci) Gur + ^ {'i^oj^ce + ^vd^ve^ci) Gae 
I ^ h. n 

DeVeie — —Heie ~ {'i'^uj^ce + ^vd^ue^ci) Gur + {XuiLO,^e ~ i^vd^ce^ci) Gne 

H — '-^^eekzGiiz, 

Velz = —Gelz H ^^—Giiz, (41) 

where 



r)oj — HiHed — ^udl^ei) 

Ser — ^ue {ViOr — "^eOr) + i^ce {ViOe — VeOe) , 
See = iOue {Vioe - Veoe) - i^^ce (^^iOr " f^eOr) • 



7. Perturbed electric current: General expressions 



In this section, we find the perturbed electric current 

, Iji'jo y -r — 

3 



Using Eqs. (38)- (41), we can write the components of this current in the form 



3lr = 



3ie = 



■ ■ rr 

rriiiA 



Wrr + Wzrr 
Wer + W,6r 



Eir + 



Eir + 



Qi 



iriiA 
Qi 



Wre + W,re 
Wee + W,ee 



Eie + WzrzEiz, 
Eie + WzezEiz, 



jiz = W.^rEir + W^^eEie + W^^^Ei^. 



The expressions for Wrr,e, Wer,e, Wzrr,e,z, Wzer,0,z, and Wzzr,e,z in the general form are given 
in the Appendix F. 



8. Dispersion relation in the general form 

Prom Maxwell's equations (4) and (5) we have, 

nlEir,e = — jir,e, = ju- 
to 



(42) 



Using the components of the perturbed electric current given in § 7, we will find from the 
system (42) the dispersion relation 

{pl^zz — ^rr^zz + ^rz^zr) {jT^z^zz — ^OO^zz + ^Oz^zej—^—^re^zz + ^rz^zo) {—^dr^zz + ^zr^Oz) — 0. 

(43) 



Here 



^Or — 



Qi 



rriiA 
Qi 



T't^l^y I "^^^ ZTT 

Wer + W,er 



1 ^rO — 



Ani 
Am 



Qi 



Qi 
rriiA 



Wre + Wzre 
Wee + Wzee 



, Ej-z ^' zrz -I 

UJ 

, ^ez — — vVzez, 

UJ 



Am Airi 

Wzzr:^z0 = WzzOt^zz = Wzzz- 

UJ U UJ 
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9. Dispersion relation in the case of the strong collisional coupUng of species 

We will further consider the case of the strong collisional coupling of neutrals with ions 
and dust grains and dust grains with neutrals. We also suppose that the collisional frequen- 
cies of electrons and ions with neutrals are much larger than the perturbation frequency. 
Thus, the case under consideration is the following: 

> ^e, > ^i, > ^d. min {i/O., I/O J > (44) 



where 

, .0 

,,0 _ '^in'^nd ,,0 _ '^dn'^ni 



The collisional frequencies Vj [Eq. (36)] are 




,,0 ,, _ ,,0 • , '^in'^ni _ • , '^dn'^nd 
^en^ ^i-^i - ^^n ^d - ~ *<^n 55—- 



Here we keep the corrections proportional to u because the main terms are cancelled when 
calculating some necessary expressions in the case (44). For example, the value A^^ has in 
this case the form 

X^^u (i/^ + + = uja, 

where i/° = t^^.t^l/i^". 

Let us find Wrrfi and Werfi under conditions (44). We will consider, as in the steady 
state, that the electrons are magnetized, uol^ 3> v^^- This condition is generally satisfied. 
Using the condition of quasineutrality QjUj^ = and carrying out the calculations, we 
obtain the following expressions for Wrr,e and Wgrfi. 

Wrr = Wee = -iqeUeoUJa, Wre = -Wer = iqefieooj a— - Ed] , (45) 



where 



'^ci 02 
= ^d^ I'd 9: 

_ Qdndo 02 f 

^d — ^ci^d J- 

qerieO 



The values g and / are equal to 



,0 '^nd 



^en ( + ) ' / i 7/0 ) ■ 
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The values Srfi,zz under conditions (44) and » are 



Srz — ^1^,7 



— {Vm - Veoe) - {ViOr - VeOr) 



QdlT-dO s 



— {ViQr - VeQr) " {vme - ^eOe) 



QdlT-dO ^ 

- J^dJ't^i09, 

qerieQ 



—UJ 



(46) 



Using the expressions for eij, i,j — r,9,z, given in § 8 and taking into account Eqs. 
(45), (46), and the expressions for Wzrr,e,z, Wzer,e,z, and Wzzr,e,z given in the Appendix F, 
we find the following relations containing in the dispersion relation (43): 



^rr^zz - ^rz^zr = T^a [ihzz - 1) + ^0^^ {hzVeOr " Kz) {ViOr " ^^eOr) , 



-EorSzz + ^ez^zr = "T/^ ( O"— " ^d ) {i^zz - 1) + ^0^^ {hz " bzzVeOe) (^'iOr " ^'eOr) , 



eeeSzz - eez^ze = rfia {ihzz - 1) + ^0^^ (hzVeoe - hz) {vm - Veoe) , 



- ^rO^zz + ^rz^ze = ^(5 ( CT— - Ed ) (ifc^z " 1) + «/?^^ (&r-z " ^^zt^eOr) ("yioe " -"60^) , (47) 



where 



2 2 
T= ^ — ^. 



Substituting the expressions (47) and Szz in Eq. (43), we obtain the following dispersion 
relation: 



Co", 



pe 



^ '.4„2,..2 , ,3 ^pe 



,2 



\ ^ 



CTl^dl 



kiwi 



(48) 



to". 



-lUJ- 



pe 



^dl 



— :-ktc^Wn = 0, 



ce 
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where 





a Uci 
-Q 9, 






qdUdo J. 

QerieO 




wl = 


{Vior - Veorf + {Viot 


9 — Veoe 



In Eq. (48) the relationship 



ib,, - 1 = — ^ 



a \ OJce 



has been used. Below, we will find the unstable solutions of Eq. (48), the growth rate of 
which is proportional to wq. 



10. Solutions of the dispersion relation 



To solve the dispersion relation (48), we will consider the cases in which the unstable 
perturbations have the long and short wavelengths. 



10.1 Long wavelength instabilities 



We neglect the terms proportional to in Eq. (48) under condition 



oj^^^ai » ky. (49) 



Then the dispersion relation takes the form 

{<^l + 4i) - ^li^Wo + iww«7 — ^^^'^^ X kiwi = 0. (50) 



For frequencies in the range 



(51) 



-20- 



the unstable solution of Eq. (50) is the following: 



1/2 



UJ 



2 



(52) 



In the case opposite to (51) the instability is caused by the presence of dust grains: 

1/3 

C^l^dl 



ujciklwl 



TT .2p7r 
exp(-z- + z— 



where p — 0, 1, 2. Note, that the sign of Sdi depends on the sign of qd 



(53) 



10.2 Short wavelength instabilities 
Now consider the case in which 



^'-^{al + el,y^'<^klc\ (54) 
Then the equation (48) will be the following: 

2 ^pe^« 1 Wl UJ^e Sdl 

For frequencies in the range 



ujci > oJSdi (56) 



we find the following solution of Eq. (55): 



, ^ ^ ^, (57) 

We see that when cti > £(ii(i/g„/a;ce) or e^i > {qd < 0) there is the streaming instability. 
In the case opposite to (56), 

LUci < UJEdl, (58) 

we obtain 



LU = —l-^—. r ^. 59 



This solution describes the instability when e^i > or 841(1'^^/ cuce) > Ci- 
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11. Discussion 

In the equilibrium and perturbation states we have not taken into account the colhsions 
between ions and dust grains, supposing that z/j„ ^ and Udn ^ i^di- These inequahties can 
be written in the form < av >i„> {miPd/nidPn) < (^v >idand < av >dn^ {Pi/Pn) < >id, 
where Pi^d,n — i^i,d,n'^i,d,n (^i ^ ^n)- Such conditions can be satisfied, if the density of 
neutrals is sufficiently high. 

Now we estimate the stationary velocities of species in the case ^ I'ni and Udn ^ 
^ni-i ^nd- These conditions can be satisfied in the weakly ionized protostellar and protoplan- 
etary disks. Then, neglecting the contribution of the radial electric field under condition 
Eor/Bo < {pn/ Pi){^K/oJci){vK/c), WB obtain from Eq. (25) 

Vnoe - ——'"K, (60) 

where vk and fix are the Keplerian velocity and frequency, correspondingly. Using Eq. (60), 
we can write the condition of the strong coUisional coupling of neutrals with ions (20) in the 
form 

"^ci ,,2 r>2 



Ki » ^K- (61) 



We see from Eqs. (60) and (61) that f°Qg ^ vk- As long as v^qq ~ f°Qg [see Eq. (27)] 
and v'^QQ ~ "^noei ^^e species rotate in the magnetized regions of the disk with the velocity 
much smaller than the Keplerian velocity. This result agrees with observations (Donati et al. 
2005). Prom Eqs. (26) and (28) we obtain the radial velocities of neutrals and dust grains 
under conditions at hand 

vlr - - -^VK- (62) 

Comparing Eqs. (60) and (62), we have v^^^g ~ ~{'^in/^ci)'^ndOr- Thus, for the magnetized 
(unmagnetized) ions, ^ {^)^in, the radial velocity of neutrals and dust grains is larger 
(smaller) than their azimuthal velocity. We see from Eq. (62) that the radial velocity of 
neutrals and dust grains is directed inward of the disk. 

It is followed from Eq. (29) that the azimuthal electron velocity v^og <S v^^q and the 
radial electron velocity VgOr — {^en/^ce)''Jnoe- 

The conditions (21) and (23) of the strong coUisional couphng of dust grains with neu- 
trals can be written in the form 



i+^fi+4 



(63) 
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Let us consider some specific parameters for the protoplanetary disk with the solar mass 
central star at r = 1 AU, where Bq ~ 0.1 G and VLk ^ '2 x 10"'^ s"^ (e.g., Desch 2004). 
We take rrii = SOrUp and m„ = 2.33mp, nip is the proton mass. Then the ion and electron 
cyclotron frequencies will be equal to ~ 32 s^^ [qi = —q^ and \ujce\ ~ 1-76 x 10® s^^ 
(the sign || denotes an absolute value). The rate coefficients for momentum transfer by 
elastic scattering of ions and electrons with neutrals are < av >in= 1-9 x 10^^ cm'^ s^^ and 
< au >en= 4.5 X 10'9(T/30 K)i/2 cm=^ g"^ (Draine et al. 1983). We take T = 300 K. Then 
the condition for neglecting collisions of neutrals with electrons has the form o;^/ (<^ci+^i>i) ^ 
1.89 X 10~^{neo/nio) (see the end of § 3). We will consider the case in which uf^ » and 
n-eo/riio ~ 10~^ (Desch 2004). Then the value must satisfy the condition < 7.36 x 10^ 
s"-*". This inequality is satisfied for the neutral mass density p„ < 2.1 x 10~^° g cm^"^ or 
n„o < 5.38 X 10^^ cm^^. We will take the ionization degree as riio = 10~^n„o (Desch 2004). 
If we set pn = 10^^" g cm~^, then we obtain = 3.5 x 10^ s~^, = 4.5 x 10~^ s~^, 

I/en = 3.65 X 10^ S-\ Une = 8.53 X 10"^° S'^ . 

The azimuthal velocity of species using the parameters given above is equal to t',°,i „o(9 — 
0.49fi^, where vk — 30 km s~^ [see Eq. (60)]. The orbiting frequency of dust grains Qd = 
0A9Qk = 0.98 X 10^^ s~^. The dust grains will be unmagnetized, ujcd ^ ^d, if their mass 
rud ^ 8 X 10^^^ g ~ 4.8 x lO^mp (for = ±qe)- The dust grains, for example, with 
density of the material grains are made 0"^ ~ 3 g cm~^ and with radius > 8.6 x 10~^ /im 
satisfy the condition of unmagnetization. At = 0.1 /im, the mass of the grain is equal to 
nid ^ 1.26 x 10~^^ g. The coUisional frequency of dust grains with neutrals has the form 
i/dn — ^■7pn''^d''^Tn/'m'd (e.g., Wardlc and Ng 1999). Using parameters given above, we obtain 
Udn — 0.53 s~^. For the solar abundance value pd ^ 10~^ p„ we have Und — 0.53 x 10~^ s~^. 

The radial velocity of neutrals and dust grains is equal to |i'd,nOr| — 133 m [see Eq. 
(62)]. This velocity is directed inward. The conditions (61) and (63) are satisfied. 

Above, we have considered only one set of disk parameters. Any other parameters can 
be treated also. 

In the case (44) of the strong collisional coupling between ions, neutrals, and dust grains 
the perturbed velocity of dust grains is of the order of the perturbed velocity of ions, v^i ~ Vji 
[see Eqs. (39) and (40)]. Thus, the dust grains can participate in the fast electromagnetic 
perturbations generated by the electrons and ions via collisions with neutrals and acquire 
the large velocities. 

Let us now consider the obtained unstable solutions. For conditions used in this section 
the growth rate 7 =Im cu of the solution (52) can be estimated in the case ai > Edi as 



7 ~ i^zr) fix- 



(64) 
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This growth rate is much larger than the Keplerian frequency as long as k^r » 1. The wave 
number must satisfy the condition [see the inequality (49)] 



k^r 



\qe\neO Pn v\ 

qiriio Pi ' 



where cai — {oJcil^pi)c is the ion Alfven velocity. For the parameters given above we have 
UiQ = 2.56 X 10^ cm~^. In this case oUpi ~ 3.84 x 10^ s^^ and, accordingly, c^i — 250 km 
s~^. So far as 7 <^ z/„j, we have the second condition: (A;^r)^^^ <^ Uni/Q.x- It is interesting 
to note that the same growth rate also exists in the collisionless regime (Nekrasov 2007). 
The growth rate (53) is larger than the growth rate defined by Eq. (64). This instability is 
possible when the density of dust grains is sufficiently large: 

^ QdndO l^ni /, \-l/2 

— - r~\ — ^ cT ■ 

"-^ci l^ePeO "X 

The short wavelength perturbations considered in section (10.2) have the following 
growth rates for conditions given above. An estimation for the growth rate (57) is 

^^/k^M'^V^^^. (65) 

V QifT'iO Pi J CAi 

This growth rate is considerably larger than the Keplerian frequency in weakly ionized disks. 
For the parameters given above we have 7 ~ 2.1 x 10~^ s~^. The condition (54) at ai > Sdi 
is the following: 

k.r > 



\qe\neO Pn Vr 
qiUiQ Pi (?Ai ' 



The growth rate (59) is larger than (65) and has the form (g^ < 0) 

\qd\ndoPn 

7 

qiTlio Pi Vni C% 

For this instability the following condition must be satisfied [see the inequality (58)] : qdndo/qeneo ^ 

J^ni/l > 1- 

Thus, in the coUisional accretion disks there are possible the streaming instabihties 
involving the dynamics of neutrals and dust grains with growth rates much larger that the 
Keplerian frequency. 



In the present paper, we did not take into account the ionization and recombination 
processes in the continuity and momentum equations which, in general, can be included 
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(see, e.g., Pinto et al. 2008; Li et al. 2008). The main sources of ionization of accretion 
disks are Galactic cosmic rays and X-rays from the corona of the central star. These and 
other chemical processes can play an important physical role in the evolution of the disk. For 
example, the dead zones of the protostcllar and protoplanetary disks may be enlivened due 
to turbulent transport of metallic ions, which are charged due to interaction with ionized 
gas in the surface layers, from the surface layers into these regions (Ilgner and Nelson 2008). 

The ionization/recombination processes in the stationary continuity equations for species 
determine the ionization degree of medium. For the weakly ionized gases with large colli- 
sional frequencies, as it is the case for us, the contribution from the ionization source terms 
to the momentum equations can be ignored (e.g., Li et al. 2008). The chemical processes 
do not influence the dynamics of medium, for example, in the case of the ideal magneto- 
hydrodynamics. In our case of multicomponent disks with strong coUisional coupling of 
species, when the dynamics of each species is considered separately, the contribution of the 
ionization source term to the perturbed continuity equation can be neglected under condi- 
tion uo > $,{nno/nio), where ^ is the ionization rate per H atom. The typical growth rates 
of the streaming instabilities considered in the present paper is 7 ~ 10^^ s~^. If we take 
nno/f^io ~ 10^, then we can neglect the ionization for ^ < 10^^'^ s^^. Note, for example, that 
in the paper (Ilgner and Nelson 2008) the range of ^ ~ 10^^^ — 10^^^ s^^ is considered. For 
cosmic rays one obtains ^ ~ 3 x 10~^^ s~^ (e.g., Li et al. 2008). Thus, for not too high ion- 
ization rates and fast instabilities studied in the present paper, the ionization/recombination 
processes can be ignored. 



12. Conclusion 

In the present paper, the electromagnetic streaming instabilities of multicomponent 
coUisional accretion disks have been studied. We have explored regions of the disk where the 
medium is sufficiently ionized and there is strong coUisional coupling of neutrals with ions 
and dust grains simultaneously. We have included the effect of perturbation of coUisional 
frequencies due to density perturbations of charged species and neutrals. This effect emerges 
when there are different background velocities of species, as occurs in accretion disks. 

We have investigated in detail the steady state and found the azimuthal and radial 
velocities of species. 

The general solutions for the perturbed velocities of species with coUisional and thermal 
effects have been obtained. We have shown that dust grains can be involved in the fast 
electromagnetic perturbations induced by ions and electrons, through strong collisions with 
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neutrals which have also strong coUisional couphng with ions. The dust grains have been 
found to acquire the perturbed velocity of ions. This effect is important for their coUisional 
coagulation and sticking. 

We have derived the dispersion relation for the vertical perturbations and found the 
unstable solutions due to different background velocities of electrons and ions. It has been 
shown that the growth rates of these streaming instabilities can be much larger than the 
Keplerian frequency. 

Electromagnetic streaming instabilities, with induced dynamics of neutrals can be a 
source of turbulence in sufficiently ionized regions of coUisional accretion disks. 
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APPENDIX A 



Solutions for v„ife and Unik 

The Fourier components v^i^ and rinik are the following (for simplicity, the index k is 
omitted here and below): 

~\~iRnlz ij^vn^rzn 



DnVnie = iR 



■nlr 



+ iRnW^vn (l — Srrn — 5||n) 



~\~iRnlz I ^un^Ozn ^^^^ ' ' 



+iRnl 



20. 



n^rzn) 



^2 / ^2 

OJ^.n (1 - (^±n) ^ - i 2Q„ - — 



" 1 X 

0- 



rOn 



— = iRnir (k - It^t^ — ke) + iRnW (h + i^^k^ + iRniz (k^ - ^k^ 



-26- 



The following notations are introduced here: 



^un = <^n + i^n^ LO^^LO- kv„o, 



= «n = {2nn/r)d{r'nn)/dr, ke = ^, 

S±n — ^rrn + ^eOni ^\\n — ^zzni — + ^\\ni 



where 1,171 = r, 9, z. 



APPENDIX B 



Solutions for Vji^ and rijik 

The Fourier components Vjik and rijik are: 

DjVjir = iFjirUJuj (l - 5eej - + iFjie [^^vj^rOj + i^^ji (l - 

+iFjiz {LOi,jSoz±j — ii^j23rz±j) , 



F>jVjiz 



iFjir {^uj^rzWj - i^325ez\\j) + iFjio {uj„j5ez\\j + «<^jl(^rz||j) 



-iFjiz 



U^j (1 - 5_Li) - ^^^^^^ - i {i^jl - ^jl) Sr0j 



—Dj-^ = iFjir ( kr - ^^^e j + iFjie ike + i-^K ) + iFju ( A:^ 



U!jlLJj2 



The following notations are introduced here: 

^3 = ^Ij (1 - - ho) - ^n^j-i (1 - ^iii) - i^i^j (^n - ^js) 5r0j, 

Fjir — Gjir + Qjlr — "^-E'lr ( 1 ~ 1^-6^—] + -^EieUr^^ + Qjlr-, 

rrij \ c / rrij c 

Fjie — Gjio + Qjie — —Eio ( 1 — n^^^ ) + —Eirng-^ + Qjw, 

m.- v c / m.- c 



9 kz 

^3 
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Fjiz — Gjiz + Qjiz — -^Eiz— + -^EiVjo h Qjiz- 

ruj u rrij c 

Ujjj = Uj + ivjiOJj = oj — kvjo, 

Q. = !^,fi:2 = {2nj/r)d{r^n,)/dr, 
kkmvlj hKvlj ^ hkzvlj .... K'^lj 

Olmj - —— , ()lz:±j - — — , dizllj - — — , d±j - drrj + Ogoj, d\\j ~ 

(jUjUJi/j UJjUJuj LujLUi/j LujUJiij 

where 1,171 — r, 9 and 

rrij "-^ ruj 

Here T±Qj and Tj|oj are the unperturbed temperatures of species j across (±) and along 
the magnetic field Bq. 



APPENDIX C 



The system of equations for perturbations of ions and dust grains 

Solutions given in the Appendix B at taking into account collisions of neutrals with ions 
and dust grains simultaneously and neglecting the thermal effects, result in the following 
system of equations for perturbations of ions and dust grains: 

DiViir + ttdrVdlr + O'dOl'^^die ~ ^dr = —H^r, 

UdO 

DiViig — ad02Vdlr + CldrVdie — bdo-^ = —Hilo, 

UdO 

^li^ilz + CtdrVdlz = iGnzUJui, 

{krttdr - keade2) Vdir + (kettdr + kradei) Vdw + ojiDi— - kb^— = -kHa, (CI) 

Tiio UdO 
DdVdlr + dirVilr + 0,i9lViW — hr = —H^ir, 

n-io 

DdVdie — (J'ie2Viir + dirVae — he— — —Hdw, 

^td^dlz + dirViiz — iGdlz^vdi 

{krttir - keai02) vnr + (keQir + Kaiei) vag + UdD^— - kb^— = -kHa- (C2) 

ndo riio 
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Here Di^a — 00'^^ ^ — <^i,di(^i,d2- The following notations are introduced above: 

Hilr — —lUuiGiir + UiiGiie, Hug = —UJi2Giir — iuJuiGne, 
Hdlr — —i^udGdlr + ^^dlGdW, Hdl0 = —ijJd2Gdlr — i^vdGdW- 

.Di Di .Dd Dd 

Hilz — —i Gilz + adr—2-Vdlz, Hdlz — —i Gdlz + O'ir—^Vuz- 

air = oJt.dfJ'i, diei = i^dilJ^i-, aie2 = i^d2fJ'i, hr = -ciir-Wr - aieiwe, he = -airWe + aie2Wr, 

ddr — ^vilJ^d^ O^dei = 'i'^illJ-d, O.d02 = 'i'^i2lJ-d, bdr = ddr^r + O-dOlWe, bde = CldrWe — ad02Wr^ 



where w = Vjo — Vdo and 





'^dn'^ni .. '^in'^nd 



IJ-i — , l^d 



APPENDIX D 

Solutions of the system of equations given in Appendix C in the case kz 7^ 0, 

kr — kg — 

The exact solutions of the system of equations given in the Appendix C for Vnr^ in the 
case kz ^ 0, kr — kg — have the form 

AViir = — {uJ^d^w + 'i'^^vi^^dl^d2) Gur — {DdUJil + fJ^ifJ^d^^dl) Gug 
+^d (A — iUiiUd2) Gdlr + IJ'd {^ud^il + ^ui^dl) G did 

-fiifidi^'^Dl^ [wr {uJyiUJyd - l^il^d + uJiiUd2) + iwg (Wivdi^a + uJuiUJdi)] kzHiiz 
- HdU^'^D'^ [wr {DdUJ^i - HiUdi^vd) + iwg {DdUJii + kzHdiz, (DI) 

AViig — {DdUJi2 + fJ'ifJ'd^d2) Giir — {oJjyd^u) + i^m^dl^d2) Gug 
—fJ'd {^vd^i2 + ^^m^^d2) Gdlr + jJ-d {^uj — ^<^i2<^dl) Gdig 

-IJ,ilJ,dUj''^D~^ [-iWr {uJud^i2 + (^m(^d2) + Wg {uj„iU„d " fJ'ifJ'd + <^i2^^dl)] kzH^z 
- HdUJ'^D-^^ [-iWr {DdUJi2 + IJ'il^di^d2) + Wg {DdUJ^i - HiUd^^ud)] kzHdlz: (D2) 

where the value A^, is equal to 
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The value A (the determinant of the system is equal to Aco'^DfDl) has the form 

A = -A^ - ul^U!iiU!i2 - (^li(^dl^d2 + i^ili^i2i^dl^d2 (D3) 
-^ii|Xd {i^i2i^dl + i^i\i^d2) ■ 

The solutions for v^irfi ai'c obtained from Eqs. (D.l) and (D.2) under substitution the 
index i by the index d, and vice versa. 

The solutions for the velocities and Vdiz are the following: 

A^filz = i^udGilz — IJ'dGdlz, 

^u}Vdiz — ^uiGdiz — l^iGiiz- (D4) 



APPENDIX E 



General solutions for the perturbed electron velocity in the case kz ^ 0, kr = ko — 

Prom the solutions given in the Appendix B, we obtain the following expressions for the 
perturbed velocities of the cold electrons: 

DeVelr = — -f^elr — /J'ei^^ueVilr — i/J'ei^^elViie — IJ-ed^ve^dlr — 'if^ed^^elVdW , 
DeVeie — —Heie + ilJ,eiC0e2Vilr — jJ'ei^veVm + ifJ'ed^e2Vdlr — jJ'ed^ueVdlO, 

i 1 

Velz = Gelz H 7— [{l^il^ed — ^udl^ei) Guz + {H'dl^ei — ^vil^ed) Gdlz] , 



where 



, .0 



APPENDIX F 

Values Wrr,e, Wer,e, W,rr,e,z, W,er,e,z, and W,,r,e,z in the general form 
Values Wrr,e and Wor^e in the general form are the following: 

Wrr = Wee = A,^ {qinioUJi^d - qdUdOl^i) + qeUeO^arl + ge^eO^^^e J^''^ A, 

2 

Wre = -Wer = i^vdi^ci {qiriioUJud - qdndO/J'i) + ^eJ^eOTf Qir-2 " ^e^^eO TT"^ A , 
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where 

The electron terms proportional to ari,2 ai'e connected with the influence of the neutral-ion 
and neutral-dust collisions on the collisions of electrons with neutrals. 

Values Wzrr,e,z-i ^z0r,e,z, and Wzzr,e,z in the general form are the following: 

2 2 
TX Tl Tl 

Wzrr Pl-^ {ibrzViOr - t^eOr) ' ^zrO = /^l"^ {iKz^iW " VeOrVeOe) , ^ zrz = i'^^rz " t^eOr) , 

2 2 
71j Tl ft 

Wz9r = fh-^ (ibozViOr " VeOrVeoo) , W,ee = /^i {ibdzViod - vIqo) , W,g^ ^ ~ ^^^^'^ ' 

7X Tl 
Wzzr = {ibzzViOr " VeOr) , VF';^^^ = {ibzzViOO " ^^eO^) , W^zzz = A (^^'^^ " 1) ■ 



Here: 



where 



br,e,zz 



^ve^rfifZZ ^ci 



1 ^lJ 



sez = {qiniQUJud - Qd'ridolJ'i) vm + {^^i^e^ee + DeVeoe) , 

Szz = {QiriiQUiud - Qd'ridolJ'i) + —■ 
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